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We analytically and numerically reveal that Majorana bound states inside the vortex core are 
distinctly spin-polarized in an odd-parity topological superconductivity classified as "pseudo-scalar" 
type of the gap function in massively Dirac-Bogoliubov-de Gennes (BdG) Hamiltonian constrained 
by strong spin-orbit coupling. This result is universal for "Dirac superconductivity" whose rotational 
degree of freedom is characterized by total angular momentum J = S + L and in marked contrast to 
spin-degenerate core bound states given by the conventional BdG Hamiltonian. The spin-polarized 
vortex core can be easily detected by spin-sensitive probes such as neutron scattering and other 
measurements well above the first critical magnetic field H c \. 
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The discovery of topological superconductor opened 
a new research avenue on superconducting states. The 
topologically-protected nature together with U(l) bro- 
ken symmetry produces Majorana fermions as gapless 
zero-energy quasi-particles at surface edges, while the su- 
perconducting gap opens in the bulk body. The emerged 
Majorana fcrmion is a counterintuitive particle whose an- 
nihilation and creation operators are identical. Such a 
unique particle has a promising role in topological quan- 
tum computing utilizing non-Abelian anyon statics. This 
fascinating novelty has highly stimulated many theorists 
and experimentalists to intensively study the topological 
superconductivity (ll-fll| . However, its research history is 
not so long as the topological insulator, and most of rich 
physics still remain elusive. 

Very recently, Sasaki et al. found that 3-dimensional 
topological insulators, Bi2Se3 and SnTe, turn into super- 
conductors with carrier doping and moreover their super- 
conducting gaps are not conventional since zero-bias con- 
ductance peaks (ZBCP's) can be detected by the point 
contact spectroscopy [lil fl3j ) . The ZBCP is known to 
be observable in not only unconventional non-s-wave su- 
perconductors but also topological superconductors [l^- 
17, [3]. A typical topological example is chiral p- 
wave superconductor such as Sr2Ru04, in which gap- 
less quasi-particles assigned as Majorana fermions in- 
duce ZBCP'sJUil]. Accordingly, ZBCP's observed in 
Cu x Bi 2 Se 3 (T c ~ 0.3K) and Sni-Jn^Te (T c ~ 1.2K) 
can be regarded to be originated from their non-trivial 
topology. 

The Majorana fermion has been suggested to ap- 
pear at not only surface edge but also vortex core 
in topological superconductors. Its emergence inside 
the vortex core has been numerically confirmed by us- 
ing Bogoliubov-de Gennes (BdG) formalism for p-wave 
triplet superconductor [iJ? ]. The vortex is basically a 
mobile object, indicating that Majorana fermion is also 
movable. The character is hopeful for the topological 
computing, since vortex manipulation techniques have 
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FIG. 1: (Color) Spin-polarized Majorana bound states. The 
red (blue) region denotes dominance of up-spin (down-spin) 
component. 



been rapidly developed in the last decade. In this paper, 
we add a new insight on Majorana bound fermion inside 
the vortex core. We analytically and numerically reveal 
that the Majorana state is spin-polarized in topological 
superconductors originated from strong spin-orbit cou- 
pling. Such a case is not particular but rather universal, 
since the spin degeneracy supposed to be kept in the con- 
ventional BdG calculations is generally broken when spin 
and orbital angular momenta are coupled. 

A minimum model of topological superconductivity 
under strong spin-orbit coupling demands two orbital de- 
grees of freedom in addition to two spin ones [J, Il2l. |25|. 
Then, the starting Hamiltonian corresponds to massive 
Dirac-type one including off-diagonal gap function, in 
which the Lorentz transformation invariance together 
with anti-commutation between two fermions gives a 
mathematical restriction on possible gap functions. As 
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described in Ref. 26[ , the allowed gap function is charac- 
terized as a scalar, a pseudo-scalar, or polar-vectors when 
considering on-site Cooper pairing. The polar- vector case 
always produces a specific direction around which the 
rotational isotropy is broken in superconducting state. 
Consequently, the restriction allows six kinds of pairings, 
which arc classified as two even-parity being not topologi- 
cal and four odd parity being topological. In CuzE^Sea, 
Sasaki et ai, theoretically checked which type of these 
pairings successfully reproduces observed ZBCP's and 
suggested that the odd-parity spin-triplet pairing is most 
likely. Another superconductor, Sni-^In^Te can be also 
examined by the same scheme. In this paper, we study 
vortex core bound-states in an odd-parity superconduc- 
tivity whose gap function is characterized as "pseudo- 
scalar" keeping the rotational isotropy. In the case, 
the zero-energy Majorana bound states emerge at vor- 
tex cores as well as surface edges, and the core bound 
Majorana fermion is distinctly spin-polarized. In other 
topological cases, the spin-polarization becomes obscure 
since the rotational isotropy is broken around the vor- 
tex center and the polarized core states are mixed each 
other. We analytically and numerically demonstrate the 
result and discuss the theoretical mechanism as well as 
its experimental detection possibility. 

An effective theory for the topological superconductor 
intrinsic to spin-orbit coupling is given by the massive 
Dirac type of BdG Hamiltonian together with Nambu 
representation for the two orbital degrees of freedom as 



H= dr(${r) Mr) ) 



H-{r) A~(r) \ / ip(r) 



A+(r) H+(r) J \ Mr) 



(1) 



where 



#±(r) = Af(-iV) + ^HW ± PoHV) 7 . 



(2) 



Here, 7* is 4 x 4 Dirac gamma matrix which satisfies the 



relations 7" = a z ® 1, 7 



i=l,2,3 



iciy ® Si and 7 



with 2x2 Pauli matrices &i in the orbital space and §i in 
the spin space, ip( r ) is the Dirac spinor, ip(r) = -0^(r)7°, 
■0 c (r) = tplj , and ip c = CM i where C(= ij 2 j°) is 
the representative matrix of charge-conjugation. The 
gap functions are defined as A~ = "f°Aij 2 and A + = 
7°(A _ )^7°, and M, P u , and e are parameters depending 
on the materials. When considering only the on-site pair- 
ing interaction, the possible gap form is reduced into six 
types of functions as seen in Table I of Ref. [26[ . These 
gap functions are classified into pseudo-scalar, scalar, and 
polar vector (four-vector) associated with the Lorentz 
transformation, 



A , A o7 5 , A n ? 



(3) 



where, Ao is a scalar constant, the Feynman slash (fx 
is defined by X^7 Mc V' an d the gap function includ- 
ing ft is characterized as unit four- vector a^. Since the 
vector type of the gap form represented by with fi- 
nite aj-components (i = 1,2,3) points a specific direc- 
tion, the rotational isotropy is broken around the di- 
rection resulting in the anisotropic quasi-particle spec- 
trum in its superconducting states, even if the normal 
state is isotropic. This induced anisotropy yields the 
angle-dependent transport conductivity. For example, 
the anisotropic thermal conductivity is a clear evidence 
of the vector type of gap given by the four- vector a M [2(| . 

From the Hamiltonian Eq. ([1]) , the correspondent BdG 
equations are given as 



j°H-{r) 7°A-(r) 
7 °A+(r) 7 °i7+(r) 



u(r) 
u c (r) 



= E 



u(r) 
u c {r) 



(4) 



Here, we note that v in the conventional cigen-state form, 
(u, v) T is given as v = ij 2 u c . 

Now, let us examine the vortex core bound states. 
Firstly, to analytically concentrate on the low-energy 
physics, we study the linearized Dirac Hamiltonian: 

ff± (r) = Mo - id xl x - id yl 2 - id zl 3 ± ^7°, (5) 

with the renormalized axes (x,y,z) — > (P\x, P^y, P3Z), 
where P v is the coefficient of the linear term of P v {— iV) 
in ^-direction, and /1 denotes a chemical potential. 
For simplicity, we do not solve the gap-equation self- 
consistently and just give the spatial distribution of 
the order parameter around a single vortex A _ (r) = 
A~/(r)e* e , where 9 denotes the polar angle around c- 
axis (vortex line), f(r) is the amplitude of the order- 
parameter (/(r = 0) = and f(r) > 0). 

At the zero energy, there is a relation expressed as 
u c (r) = i-f 2 u*(r), and the solution is expressed as u(r) = 
£(r)u N (r), where u N (r, 9, z) is a solution in the normal 
states and C(r) is a scalar real function. In the case of the 
pseudo-scalar type of gap function A2 = 7 2 7°, there are 
two bound-state solutions at the zero-energy expressed 
as 



-K(r) 



( 



Ui{r) 



M J {f) 





-K{r) 



\ ie 
( 



-M Ji(f) / 
\ 



-iy/fJ. - 300 Jo (f) 

V J 



(6) 



(7) 



where f = r\/ fi 2 — M5 2 , K(r) = J Q r \f(r')\dr' and J„(r) is 
the Bessel function of the first kind (See the supplemental 
materialj^l. Here, X± is determined by 

/•OO 

A± = 4tt / drre~ 2K ^ [(/x ± M ) J Q (f) 2 + (fi T M ) Ji(f) 2 
Jo 

(8) 
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FIG. 2: (Color online) The radial dependence of the spin- 
resolved LDOS's with the gap-amplitudes (a) /i = 0.4eV and 
A = 0.3eV, and (b) fi = 0.5eV and A = O.OleV. 



We note that these two solutions are localized when 
A/ 2 < fi 2 + |Ao| 2 [13|. Totally, we can obtain the so- 
lution at finite energy by a perturbation in terms of 
k z . Assuming the solution ip(r) T — c^(u^(r),u c ^(r)) T + 
c^{u^{r), u c i(r)) T , the dispersion relations are simply 
given by 



E = ±vk z 



(9) 



where the coefficients (ef,cj.) = (l,±£)/v2 are not de- 
pendent on k z , and 



V = 27r vV_M 2 j X drre -2K(r) (j()(f) 2 _ Ji( - )2) _ 

y/X + X- JO 

(10) 

The spin-resolved local densities of states (LDOS's) 
n t (E = vk z ,r) = l^iWI 2 + \Mr)\ 2 and n^E = 
vk z ,r) = \ip2(r)\ 2 + \ip4(r)\ 2 are expressed as 

n^(E = vk z ,r) = I — — + — — J M r ) e ( , 

(11) 

/ , \ ( M - M o , M + T ,-,2 -2KM 

ni(E = vk z ,r) = I — — + — ^ — \ Jiif) e { >. 

(12) 

Here, ')/'i( r ) is the i-th component of the array solution 
form i/j(r). We show r-dependence of n-|- and n± in two 
parameter sets with the assumption of the spatial form 
of the order-parameter f(r) = r/\/r 2 + 1. As seen in 
Fig. [2fa), rif in /x = 0.4 eV and Ao = 0.3 eV shows the 
peak structure at the vortex center (r = 0) in contrast to 
the zero density in because of Jo(0) ^ and Ji(0) = 
0. In /i = 0.5 eV and Ao — 0.01 eV, which are the 
same as that in the previous papers (l2l . |26| , and 
show qualitatively same similar distribution patterns as 
shown in Fig. [2jb) . The integrated magnetization is also 
calculated by 

M z (r) = J dk z (n t (k z ,r)~n t (k z ,r)) (13) 
~ m(n t (£J~ 0,r) - nf(E ~ 0,r)), (14) 



where m is a constant. Then, it is clearly found that the 
vortex core is spin-polarized as found from Fig. [2] 

Next, we numerically calculate the LDOS's with the 
use of the material parameter sets for Cu^E^Sea to com- 
pare with the analytical results. We employ L x x L y tri- 
angle lattice grid. A single vortex center is located at 
(ix,iy) — (L x /2, Ly/2). To obtain the LDOS njuj, i X] i v ), 
we use the spectral polynomial expansion scheme [3ll433l | 
with 40 fc z -points and L x = L y = 96. We use renor- 
malization factors a = 30 eV, b = — /z, a smearing factor 
X] = 1 X 10 -3 eV, and a cut-off parameter n c — 8000 
(see, Ref. (3l|). We take the gap-amplitude A = 0.3 
eV and the chemical potential fi = 0.4 eV, which are 
the same as the parameter sets of the analytical result 
in Fig. EJa). Other parameters are the same as those 
in Ref. [26J . In the presence of a vortex, one can find 
the zero-energy bound states at the vortex core, in which 
there are the energy spectra of two bound states in the 
fcz-resolved LDOS's as shown in Fig.(3ja). One shows the 
linear dispersion relation around the zero-energy, while 
the other has a flat dispersion relation at the mid-gap 
energy (E/A ~ 0.5). We confirm that the fc z -dispersive 
spectrum consists of only up-spin quasiparticles while 
the mid-gap spectrum are opposite (See, Fig. [3fb)). As 
shown in Fig. [4j the r dependence of the spin-resolved 
LDOS at the zero-energy around a vortex shows that the 
vortex core is spin-polarized, which is consistent with our 
analytical calculation shown in Fig. [Ha). 

Finally, we discuss the reason why the majorana 
fermion inside the vortex core is spin-polarized. First 
of all, we recall that the quasi-particle energy spectrum 
of the core bound state in topological superconductiv- 
ity such as triplet p-wave one is given as El oc L 
where L is the orbital angular momentum in contrast to 
El oc (L + 1/2) in non-topological ones. Here, it should 
be noted that the above relation is valid when the full gap 
opens and the cylindrical symmetry is preserved around 
the vortex center. The above relation clearly indicates 
the presence of the zero-energy Majorana bound state at 
L = in topological superconductivity. Then, we notice 
that the Majorana state is degenerated in terms of its 
spin momentum S in the conventional BdG scheme, re- 
sulting in no spin-polarization. On the other hand, when 
the spin-orbit coupling turns on, the preserved angular 
momentum L is replaced by the total one J = L + S . 
Then, a core bound state without the spin-orbit coupling 
is splitted by the spin degree of freedom into two states 
whose total angular momenta differ by 2S(— 1) and the 
mostly-bound Majorana state at the zero energy is a state 
of the splitted ones. Thus, one finds that the bound 
Majorana particle is distinctly spin-polarized. This con- 
sequence is universal and analogous to atomic orbitals 
when increasing the nuclear charge Z. The spin-polarized 
core is easily observable if vortex lattice is formed. For 
example, the neutron or muon scattering is sensitive to 
such a spin-density lattice like modulation, and NMR is 
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also a good probe. 

In conclusion, we obtained vortex core bound-states 
in a pseudo-scalar type topological superconductor using 
the model Hamiltonian of the topological superconductor 
Cu K Bi2Se3. We found that the zero-energy Majorana 
fcrmions are spin-polarized and the vortex core itself is 
magnetized. Such a drastic feature is easily observable 
in experiments. 

We thank M. Okumura for helpful discussions and 
comments. The calculations have been performed using 
the supercomputing system PRIMERGY BX900 at the 
Japan Atomic Energy Agency. 
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FIG. 3: The energy dependence of the LDOS's at the vortex 
center (i x ,i y ) — (L x /2, L y /2) in the presence of a vortex in 
the superconductor whose the gap function A2. (a) the k z - 
resolved LDOS's. (b) the total and partial LDOS's. 
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FIG. 4: The r dependence of the spin-resolved LDOS's at the 
zero-energy, r is determined by r = V3(i x — L x /2)/(2A2). 
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Supplemental Material 



SI. Solutions of the BdG equation around a vortex 

Let us solve the following BdG equations: 

H 



where the Hamiltonian is 



with 



H 



u(r) 
u c {r) 



7°#-(r) 7°A-(r) 
7 °A+(r) 7 °i7+(r) 



E 



u{r) 
u c (r) 



h(r) - n A'(r) 
A' f h(r)+fj, 



h(r) = 7°M - ifl^V - z^ 7 °7 2 - i0*7°7 3 



(1) 



(2) 



(3) 



It should be noted that one can find the solution (u , u c ) = (i7 2 u*, i7 2 u*) with the energy — when (u,u c ) is a 
solution with the energy E. The zero-energy solutions satisfy the following equations, 



(h - n)u + A'i 7 V = 0, 
A'Uj 2 u* c + (h + fi)u c = 0. 

We assume that these solutions are expressed as 

u(r) = ((r)u N (r,9,z). 

Here, u N (r, 6, z) is a solution in the normal states which satisfies 

(h-n)u N (r) =0. 

At first, we solve the equation in the normal states. The equation is rewritten as 

h 2 u N (r) = f J 2 v N (r), 
(-V 2 + M 2 )u N (r) =n 2 u N (r). 

Thus, each component of the solution is expressed as 

u m (r) = Ci e inid J ni (ar)e ik * z , 



(4) 
(5) 

(6) 
(7) 



(8) 
(9) 



(10) 



with a — \/ [J? — Mq — k 2 . In order to obtain coefficients (01,02,03,04), we solve the following equations, 



( M k z L_ \ 

M L+ -k z 

k z L_ -M 

V L+ -k z -M j 



/ Cl e m ^J ni {ar) \ 
c 2 e m ^J n2 {ar) 
c 3 e m - e J m (ar) 

\ C4 e m ^i 4 (ar) / 



Ah z z 



( Cl e in ^J ni (ar) \ 
c 2 e m * e J n2 (ar) 
c 3 e m ^ e J n3 (ar) 

\ Ci e m * e J ni (ar) J 



(ii) 



with 



d 1 d 



dr r 89 

With the use of the relations d r J n (r) — —J n+ \{r) + (n/r)J n (r) and d r J n (r) = J„_i(r) — (n/r)J n (r) , we obtain 

L ± e m0 J n (ar) = ±e« n±1 > fl (Ti£ + -)Uar) 

or r 

= ±e l( - n±1 ^iaJ n±1 {ar). 



(12) 

(13) 
(14) 
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Thus, the equations are rewritten as 

/ (M -/x)e ini9 J ni (ar) k z e m * e J n ,(ar) -e^-^aJ^-^o 

(Mo-//)e in ' 9 J n2 (ar) e l ^ +1 ) e m J„ 3+ i(ar) -k z e in * e J„ 4 (or) 

fc 2 e mif, J ni (ar) -e'^-^W^^far) (-M Q - A*)e m3e J„ 3 (ar) 

V e l ^ +1 ^taJ ni+1 (ar) -k z e m * e J„ 2 (ar) (-Af - M )e m ^J„ 4 (< 



The orbital angular momenta n% must satisfy 

n 2 = ni + 1, "3 = Tii, tl 4 = ni + 1. 
Then, in the case of k z — 0, the coefficients are written as 

ic 2 a ic\a 





(ci\ 




C2 




C3 


) 





C3 



M + ^ : 



C4 



M + /x' 



= 0. 
(15) 

(16) 
(17) 



We assume that £(r) in Eq. (|6]) is a real scalar function. By substituting Eq. (J6j> into the differential equations Eq. Q, 
we obtain 



A r ie dC(r) ic 2 a in^B j ( \ 
' e dr Mo+fJ, J ml"'; 

ie^ s ^c 2 e^+^J ni+1 (ar) 



V 



+ C(r)A (r) 



^e-^ e J ni (ar) 



V 



„ c * e -,(« 1+ i)e Jni+i(ar) 

? J ni (ar) / 



(18) 



Finally, we consider the gap function represented by s pseudo-scalar. In this case, A'(r) with a vortex is expressed 



A'(r) = /(r)e* M V 



(19) 



Here, M is a winding number and f(r) is the amplitude of the order parameters (f(r = 0) = and /(r) > 0). By 
substituting A'(r) into the above differential equations, we obtain 



; d((r) 



c 2 e m ^ J„ 1+ i(or) 



C0)/(r) 



The angular momentum of the zero energy solutions is 



^ S -n 1+ M)ej nMr) 

-c* 2 e^ n ^ 1+M ^J ni+1 {ar) 
V c* 1 e^-^ +M ^ 6 J ni (ar) J 



M - 1 



ni 



Since £(r) is a scalar real function, we obtain two solutions with coefficients 



(ci,c 2 ) = (|ci|exp 



,0), (0, |c 2 |exp 



m 

T 



)• 



Therefore, the zero-energy solutions with a winding number M are written as 

( + Aln exp [i^6] J M -i(y/u 2 - Mfrl \ 



\^ iy/n - M exp [i^!0] Jw (\/u 2 - Aigr) J 

( ° 

+ M exp jM+ifJut-Mgr) 

-Wfi - M exp [i^fl] Jm-iU/p? - MZr) 



withlf(r) =/ r |/(r')Mr'. 



V 



/ 



(20) 



(21) 



(22) 



(23) 



(24) 



